Heat conduction in the diatomic Toda lattice revisited 
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The problem of the diverging thermal conductivity in one- 
dimensional (1-D) lattices is considered. By numerical simu- 
lations, it is confirmed that the thermal conductivity of the 
diatomic Toda lattice diverges, which is opposite to what one 
has believed before. Also the diverging exponent is found to 
be almost the same as the FPU chain. It is reconfirmed that 
the diverging thermal conductivity is universal in 1-D systems 
where the total momentum preserves. 
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Heat conduction in one-dimensional lattice is rather an 
old problem. Many authors have investigated the prop- 
erty of thermal conductivity to understand what ingre- 
dients are essential to the standing of a macroscopic law, 
that is, the Fourier law; 



(j) = -k V T, 



(1) 



where k is a thermal conductivity. It is well known that 
in integrable systems such as harmonic chains or ideal 
gas, the Fourier law is not valid since no temperature 
gradient is formed Q, while various numerical simula- 
tions of nonintegrable systems show temperature gradi- 
ents. However, it is also found that thermal conductivity 
of nonintegrable systems such as the FPU chain diverge 
as N a HH where N is the degree of freedom. In other 
words, thermal conductivity becomes infinite in the ther- 
modynamic limit. 

On the other hand, finite conductivities are seen in 
some 1-D nonintegrable systems. Casati et al. invented 
so-called ding-a-ling model consisting of alternate har- 
monic oscilators and free particles, and found that the 
model has the finite conductivity g. The similar kind 
of the model which has finite conductivity is also inves- 
tigated by Prosen and Robnik ||. Most recently, Hu et 
al. found that the Frcnkcl-Kontrova model has the finite 
conductivity ||. In these models, the conductivity con- 
verges at the certain value with relatively small N which 
does not exceed 100. This convergence makes apparent 
contrast with the FPU chain where the conductivity still 
grows even at N - 5000 §,§. As the common feature of 
these systems which have finite conductivity, the exter- 
nal field is introduced to confine the movement of each 
particle. The Hamiltonian of the systems are represented 
as 



H = 



— + U(x i+ i - Xi) 
2m 



V(xi)], 



(2) 



point, one might think that the external field plays the 
key role for obtaining the finite condutivity ||. However, 
the finite condutivity is also obtained for the diatomic 
Toda lattice (DTL) m whose Hamiltonian is written as 



H 



E 



2m,j 



+ exp(.T 4+ i - Xi)], 



(3) 



where rrii denotes the mass of alternate two different par- 
ticles. The DTL has no external potential which is dif- 
ferent from Eq.(||). It has been still under cover what is 
responsible for the finite conductivity. 

Recently, Lepri et al. found that the autocorrelation 



function of the total heat current vanishes like t 



-0.6 - 



. the 



FPU chain || . This implies the divergence of the thermal 
conductivity as a result of the Green-Kubo formula; 



k = lim lim 

t^oo y^oo Vk B T 2 



dt'{J(0)J{t')), 



(4) 



where J(t) = J j(x,t)dx, and V is the volume of the 
system. 

Indeed, due to the conservation laws, long-time tails 
of the correlation functions are quite general results in 
fluids ju^. The rough explanation is as follows. Hydro- 
dynamically, the local heat current j{x,t) is expressed 
as 



j(x, t) = h(x, t)v(x, t) — k y T(x, t), 



(5) 



where V(x) is the external trapping potential. At this 



where h(x,t) and v{x,t) denote local enthalpy density 
and local velocity of the fluid, respectively [jllj] . Since 
v(x,t) appears in the first term of Eq.(||), an autocorrela- 
tion function of the total heat current (J(0) J(t)) includes 
effect of the velocity autocorrelation function (VACF) 
]l2[ . In the system where the total momentum pre- 
serves, asymptotic behavior of the VACF is propotional 
to t~ d / 2 , where d is the dimensionality of the system. 
Then (J(O)J(t)) also decays like t~ d / 2 , which implies the 
divergence of the integral of Eq. for d < 2. In the 
system where the total momentum does not preserve, the 
VACF vanishes much faster than that. For example, in 
the Lorentz gas the VACF decays like — t~ d l 2 ~ x [y_3|. In 
those systems the VACF does not cause the divergence of 
Eq. ([|) . We remark that the contribution of the second 
term of Eq. (§) to (J(O)J(t)) is r 11 / 2 - 1 . This term is 
not responsible for the diverging conductivity. 

Those explanations account for the diverging conduc- 
tivity in the FPU chain, and also the finite ones of the 
models where the total momentum does not preserve due 
to the external field, such as the ding-a-ling model. How- 
ever, the explanation does not apply to the diatomic 
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Toda lattice where the total momentum preserves. The 
fact that the DTL has a finite thermal conductivity has 
been invoking confusions. In this Rapid Communication, 
we recheck the result of Ref. 0] to find out what is really 
going on in the DTL. 

The Hamiltonian of the DTL is given by Eq. ||). 
We perform numerical simulations of the DTL in con- 
tact with two thermal reservoirs whose temperatures are 
denoted as T\ and T%. Note that the choice of models for 
thermal reservoirs is critical, since there might exist the 
temperature gaps at the extrema of the lattice connect- 
ing with the reservoirs. It makes the definition of tem- 
perature gradient ambiguous, because the system will not 
obey the assigned boundary conditions; i.e. temperatures 
of the thermal reservoirs. Since thermal conductivity is 
defined as (J)/ V T, it is important to determine \jT 
exclusively by controll parameters. The model we adopt 
here is the thermal wall type |4[|l4[]. When the particle 
collides with the wall, it reflects the particle back with a 
new momentum p at random. The probability distribu- 
tion function of p is given by 



\p\ p 2 



rnksT 2mkBT 



(0) 



The local heat flux ji(t) is defined as the energy trans- 
fer per unit time from the 1-th particle to the (l+l)-th 
particle. 



dU(xi - xi+x) 
Jl(t) = tei V " 



(7) 



The total heat current appearing in the Green-Kubo for- 
mula is 



J(t) 



N 

E 

i=i 



ji(t)a, 



(8) 



where a is the average distance between two particles. 
The average current is then defined as 



(9) 



Hereafter we fix the mass ratio of the particles to be 
0.5, that is, m,2 n — i = 2r7j 2n . The temperatures of the 
thermal reservoirs are set to be 100 and 10. Note that 
all these conditions are the same as Ref. || except for 
the reservoir model. Numerical integration is done by 
the symplectic integrator of the fourth-order |jl5|| in order 
to preserve the symplectic structure of the phase space. 
Note that the distance between two thermal walls is aN 
so that the average density is fixed regardless of the num- 
ber of particles. We set a = 1 and m 2 „ = 1 for non- 
dimensionization. 

First we check the temperature profile. We define the 
temperature of the 1-th site as the long time average of 
mivf based on the virial theorem. The result is shown in 
Fig. 1. Since no gap is seen at the extrema, temperature 



gradient \/T becomes (Tj — T^j/N. We can safely define 
the thermal conductivity as 



U)N 
T\ - T 2 



(10) 



where (J) is defined by Eq.(^). The system size depen- 
dence of the thermal conductivity is shown in Fig. 2. It 
is clearly seen that the conductivity diverges like jV - 35 . 
The exponent 0.35 is very close to the one for the FPU 
chain (0.38). It is reasonable to consider that the origin of 
the divergence is the same as the case of the FPU chain, 
i.e. the long-time tail of the Green-Kubo integrand. We 
check an autocorrelation function of the total heat cur- 
rent (J(O)J(i)), by taking a periodic boundary condition 
instead of thermal walls. The initial condition is chosen 
within the microcanonical ensemble whose temperature 
is (Tj + T 2 )/2 = 55. Fig. 3 clearly shows the long-time 
tail which is approximately propotional to £~ - 65 just like 
the FPU chain. This long-time tail is the strong evidence 
for the diverging thermal conductivity in Fig. 2, and also 
helps the unified understanding of the heat conduction in 
1-D lattices. 

However, one may think that the temperature differ- 
ence adopted here is so large that the linear response the- 
ory does not apply. To answer the suspicion, we check the 
thermal conductivity at the smaller temperature gradient 
which is closer to equilibrium, i.e. T± = 5, T 2 = 4. The 
system at this temperatures also shows the divergence of 
N - 35 and the long-time tail of <-°- 65 . 

In order to confirm the divergence in the diatomic Toda 
lattice, we also test other versions of the DTL [[16i ; 



Ep 2 
— h hardcore, 
2m l 



(11) 



H = y^[^~ + ex P(^i - x i+i) + x i+i - Xi]. (12) 



For thermal reservoirs, we use thermal wall model as be- 
fore in the diatomic hard spheres of Eq. ([ll]) . Note that 
thermal reservoir employed in the simulation of Eq. ([H 
is the Langevin type. 

miii + Qvx + = 1 - exp(xi - x 2 ), (13) 
m N v N + C"JV + &(*) = -1 + cxp(a; A r_i - x N ), (14) 

where denotes the Gaussian white noise. — 
and (&(0)&(t)) = 2Ck B T l 8(t).) We set T x = 5 and 
T 2 = 4 for the both model. System size dependences of 
the thermal conductivity of these models are shown in 
Fig. 2. They also show the divergence of TV 33 ~ TV - 37 . 

The result obtained in this Rapid Communication is 
quite opposite to the results of Jackson et al. ||. The 
keypoint is the formation of the temperature gradient. 
In Ref. ||, the temperature profile has large gaps at the 
extrema of the lattice so that the real temperature gradi- 
ent gets smaller than N/(T\ —T2). Hence it is improper 
to define the thermal conductivity as {j)N/(T\ — T 2 ) as 
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they did. Moreover, since the size of the gap may depend 
on N, system size dependence of the thermal conductiv- 
ity measured in that way is not precise. 

The existence of the gaps is due to the model of the 
heat bath. In Ref. J8|, the new momenta are randomly 
given to the end particles of the lattice. Although the 
distribution function is the same as ours, i.e. Eq. (|^), 
the new momenta are given at finite time steps which 
is determined randomly from the uniform distribution. 
When the average time interval is shorter than the relax- 
ation time of the lattice, the gap is formed. This issue 
has been partially reported in Refs. In our mod- 

els, for instance, the Langevin model represented by Eqs. 
( |l3| ) and (|l4|) yields temperature gaps when ( becomes 
large. 

In this Rapid Communication, we have confirmed that 
the thermal conductivity of the diatomic Toda lattice 
diverges as TV 35 just like the FPU chain. This diver- 
gence is gereric in the 1-D momentum preserving sys- 
tems, due to the long-time tails in the Green-Kubo inte- 
grands. Only are the systems where the total momentum 
does not preserve and the 3-D fluids expected to have the 
finite thermal condutivities in the thermodynamic limit. 

However, it is still unclear that the quantitative con- 
ditions for the existence of temperature gradients, aside 
from the choice of the heat bath model. Nonintegrability 
itself is the necessary condition. Quantitative study of 
the transport processes from the viewpoint of dynamical 
systems must be the main focus of the future problem. 
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FIG. 1. Temperature profile formed in the diatomic Toda 
lattice. The temperatures of the reservoirs are 100 and 10. 
System size N is 1000. The shape of the profile will not 
change with the increase of N. 

FIG. 2. System size dependence of the thermal conduc- 
tivity. Circles correspond to the diatomic Toda lattice with 
Ti — 100, T2 = 10. Squares denote diatomic hard spheres 
of Eq. (11). Triangles represent another version of the DTL 
written as Eq. (12). The solid line is propotional to jV ' 35 . 

FIG. 3. The autocorrelation function of the total heat cur- 
rent with the periodic boundary condition. Dashed line is 
propotional to t~ 65 . System size N is 2000. 
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